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Abstract 



Exact density profiles in the steady state of the one-dimensional fully asym- 
metric simple-exclusion process (FASEP) on a semi-infinite chain are obtained 
in the case of forward-ordered sequential dynamics by taking the thermody- 
namic limit in our recent exact results for a finite chain with open bound- 
aries. The corresponding results for sublattice parallel dynamics follow from 
the relationship obtained by Rajewsky and Schreckenberg [Physica A 245, 
139 (1997)], and for parallel dynamics from the mapping found by Evans, 
Rajewsky and Speer [J. Stat. Phys. 95, 45 (1999)]. Our analytical expres- 
sions involve Laplace-type integrals, rather than complicated combinatorial 
expressions, which makes them convenient for taking the limit of a semi- 
infinite chain, and for deriving the asymptotic behavior of the density profiles 
at large distance from its end. By comparing the asymptotic results appro- 
priate for parallel update with those published in the above cited paper by 
Evans, Rajewsky and Speer, we find complete agreement except in two cases, 
in which we correct technical errors in the final results given there. 
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I. INTRODUCTION 



We consider the current and density profiles in the steady state of the fully asymmetric 
simple-exclusion process (FASEP) with open boundaries and different discrete-time updates, 
namely, ordered sequential, sublattice parallel and parallel. We remind the reader that the 
model describes a system of particles on a chain, hopping with probability p only to empty 
nearest-neighbor sites to the right. Each of the L sites of the chain can be either empty or 
occupied by exactly one particle. Open boundary conditions mean that at each time step 
(update of the whole chain) a particle is injected with probability a at the left end of the 
chain {i = 1), and removed with probability (3 at the right end {i = L). The definition of the 
model includes the choice of the stochastic dynamics, i.e., the update scheme which specifies 
the order in which the local hopping, injection and particle removal are implemented. The 
case of random- sequential update is considered as a lattice automaton realization of the 
corresponding continuous-time process. In the general case it was solved first by using the 
recursion relations method 0, and then by means of the elegant matrix-product Ansatz 
(MPA) [Q. As it was proven later IQ, the MPA is actually not an Ansatz, but an exact 
representation of the stationary state of any one-dimensional system with random-sequential 
dynamics involving nearest-neighbor hoppings and single-site boundary terms. The method 
of MPA was next successfully applied for solving the following basic cases of true discrete- 
time dynamics: forward ordered sequential (— >), backward ordered sequential (^) p, 0, 
and sublattice-parallel (s- ||) 0], |^, which turned out to be closely related 0. Thus, the 
current has the same value in all these cases, Jj^ = = , the local densities at site 
z G {1, . . . , L} for the FASEP with forward, p'^ii), and backward, p'lii), ordered sequential 
updates are simply related to each other, pj^{i) = — J]^, and to the local density 

"(i) for the sublattice-parallel update 

pV\^) = i^ifl' (1) 

[ Pl [t), t even. 

In the cases of random-sequential, ordered sequential and sublattice parallel updates, 
the matrix-product representation involves infinite dimensional matrices which satisfy a 
particular quadratic algebra. Finally, the most difficult case of fully parallel dynamics and 
open boundaries was solved by Evans, Rajewsky, and Speer |]10|, by using site-oriented MPA 



with matrices satisfying a quartic algebra. Moreover, these authors showed that the current 
and local densities for the model with parallel update can be simply mapped onto those for 
the previously solved discrete-time updates. An independent, bond-oriented MPA solution 
for the stationary FASEP problem with parallel dynamics was found by de Gier and Nienhuis 



111 . In the case of general values of the hopping probability p, they have presented explicit 
expressions for the current and the discrete slope tiii) = Pbii + 1) — Piii) of density profile 
in all qualitatively different domains of the parameter space. 

In our recent paper |]r2[ we have derived, independently of [^, |TT|], and by using a 



different method, exact expressions for the steady state current and the local density 
p'lii), i G {1,...,L}, for the FASEP on a finite chain with forward-ordered sequential 
dynamics and open boundaries. These expressions involve integrals which at large L and 
large distance from the chain ends are of the Laplace type, hence, they are convenient for 
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asymptotic analysis. We point out that from the simple mapping [|T0| of the above quantities 



onto the current J I and local density p\{i) for the model with parallel update, 

exact new representations for j| and pi(^) follow, which do not contain complicated com- 
binatorial expressions as those obtained by Evans, Rajewsky, and Speer. 

We shall present here exact expressions for the local density profiles of the model on a 
semi-infinite chain with either left-hand (1) endpoint, 

pZ,\ii\(^^l3) = lim pL{i\a,l3), (3) 

or right-hand (r) endpoint, 

pZA^l^^f^) = /™ Pl{L~ J + l\a,(3). (4) 

In the former case the limit L ^ oo is taken at fixed i = 1,2, . . ., which labels the sites 
of the chain from its left end to the right. In the latter case the limit L ^ oo is taken at 
fixed j = 1,2, . . ., which labels the sites of the chain from its right end to the left. Note 
that for fixed hopping probability p, and any point {a,j3) on the phase diagram, see Fig. 1, 
the right-hand profile p'^j.{i\a, P) is closely related to the left-hand profile for 
the point which is symmetrically positioned with respect to the diagonal a = p. Indeed, 
by taking the L ^ oo limit in the symmetry relation for a finite chain 0, p]^{i\a,j3) = 
1 — Jj^ — pj^ {L — i + l\/3 , a) , one obtains that 

P^,i(^l"' /5) = 1 - -^oc^ - PZA^lf^^ (5) 

Therefore, it suffices to derive explicit expressions for the local densities p'^i{i\a,P) and 
p^^{j\a, P) for, say, a < /?, i.e., in the maximum-current and low-density phases. Moreover, 
by applying the L — > oo limit form of the mapping (H) with i fixed at a finite distance from 
one of the chain ends (e), 

JL = T^^ pUO = ^%^j^, (e = l,r), (6) 



to the results for p^ g(i), one obtains the corresponding expressions for p|J^ ^(z), i = 1,2, . . 
e = l,r. We emphasize that the resulting expressions for the density profiles on semi-infinite 
chains are much simpler than those on finite chains, and they again involve integrals which 
are of the Laplace type for large distances from the chain end. Hence, their asymptotic 
behavior is readily deduced. In order to make precise comparison with the asymptotic 



expansions for the right-hand local density profiles obtained in |]T0[, one has to take into 
account that there, before taking the limit L ^ oo, the sites are labeled from the right end 
of the chain to the left by n = L — i = 0, 1, 2, . . ., rather than by j = L — i + 1 = 1,2,..., 
as in our definition (^). By comparing the large j asymptotic behavior of pLrO); obtained 
via the mapping (|^), with the corresponding results of Evans, Rajewsky and Speer |10 



for pn with n = j — 1, we find complete agreement in all but two cases, in which we 
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correct technical errors in the final asymptotic expressions (9.22) and (9.28) given there. 
The expressions for the slope of the density profile derived in |ll]] for general values of the 
hopping probability p actually represent the leading-order asymptotic form of tL^i) that 
follows from our expressions for the local density Pl(^) HH in the semi-infinite chain limit 
L ^ oo and in the leading-order expansion in finite but large L — i = n ^ 1, see the results 
and comments below. For recent general reviews on the MPA, ASEP and related systems 
we refer the reader to and |T3|. 

Concerning the notation, we make the following remarks. In contrast to our previous 
work here we shall adhere to the conventional labeling of the different regions in the 
phase diagram, see Fig. 1. The mean field line shown there is given by the equation 
{1 — a){l — (3) = 1 — p, and the two other special lines are defined by 



a = ar 



1 



(7) 



The exact finite-chain results obtained in [0 are conveniently expressed in terms of the 
parametrs 



d = \Jl—p, a = d + d ^, ^ 
which will be used here too. 



p — a P — P 

ad ^ ^ f3d ^ 



(8) 



II. LOCAL DENSITY IN THE MAXIMUM-CURRENT PHASE 

The maximum- current phase occupies the region and /3c < < 1, see region 

C in Fig. 1. In terms of the variables (^ the above inequalities read: —d < ^ < 1 and 
-d <T] < 1. 



The final exact result obtained in |12| for the local particle density in the maximum- 
current phase is 



Pl(0 = ^ (1 - Jl) + ^[i^L(^) - (e - V)Z^^lZL^ 

2 2pZL 



(9) 



Here Zn = Zn{C,,ri) has the representation ^ rj) 
which involves the integral 



'd' 



^ 4(0 + ^/„(r^) 



(10) 



4(0 = - fd^^ 

TT JO 1 



a + 2 cosf^)"" sin 



(11) 



The expression for Zn{^,C,) follows by taking the limit 77 ^ ^ in Eq. ([TUp. The current 
Jii^^v) is given by 



M^,r^) = ZL-i{^,r])/ZLi^,r^). 



(12) 
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The term -^^(i) = FL{i;^,r]) in Eq. (j^) is an antisymmetric (with respect to the center of 
the chain) function of the integer coordinate i, F^ii] ^, 77) = —Fl{L — i + 1; ^, r^), defined for 
1 < i < [L /2] {\x\ denotes the integer part of x) by the equation 



.p. 



I E^/L-.-n-l(0Wl(^)■ 



(l3) 



n=0 



In order to take the hmit L ^ 00 in Eq. 
form of the Laplace- type integral (PH): 



we make use of the large-n asymptotic 



IniO = 2^/^(1 + Oin~-% (e 1) 



(a + 2) 



n+1/2 



-n 



-1/2 



[l + 0{n-% 



(14) 



By substituting the above expansion in Eq. (|TU|), we obtain that in the maximum- current 
(m.c.) phase the leading asymptotic form of Zn{^,ri) is (^ 7^ 1, 77 7^ 1) 



'dy (a + 2)-+3/2 



2v/?(l-e)2(l-r/)2 \pj 
Hence, the asymptotic form of the current reads 



3/2 



[l + 0{n-% 



(15) 



jm.c. 



[1 + 0{L-'] 



1 + v^T^ 

independently of the parameters a and f3. With the aid of Eq. (p!5D , at fixed i we obtain 



(16) 



lim 

L^oo 2p 



2(e-r7)(a + 2) 



(17) 



Next, we split the sum over n in Eq. (|T3]) into two parts, from to [L/2] — i and from 
[L/2] — i + 1 to L — 2i, and replace the cofactor in the product lL-i-n-i{0^i+n-i{v) which 
has large subscript as L ^ 00 by its asymptotic form (0). Thus, we prove the limit 



d FL{i;^,r]) 



lim 



-| 00 
2 ^ 

^ n=0 



:i-o'Wi(o^(i-^)'Wi 



(a + 2) 



(a + 2)^+" 



+1 



(18) 



Finally, by using the identity 



1 g (1 - _ G,^^ - e/i_i(0 



n=0 



(a + 2)*+"+i 



(a + 2)'' 



where 



Gn — — 



TT JO 



d0(a + 2cos0)"(l + COS0), 
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(19) 



(20) 



and collecting the above results, we obtain an exact expression for the local density profile 
of a semi-infinite chain with left-hand endpoint in the form 



d 



l + d (a + 2) 



-[G,_i-^/,_i(0]. 



(21) 



The approach of the left-hand profile for finite chains, obtained by computer simulations, to 
the limit (^) as the chain size increases is illustrated in Fig. 2. 

The asymptotic behavior for z ^ 1 readily follows from the expansions ([I^) and 



+ 2)*-V2 r 



TT 



-1/2 



-3/2 



Thus we obtain 



d 



Poo,M 



+ 



16 



-1/2 



+ 0(z 



-5/2^ 



(22) 



l + d ^{l + d) 



v^(l + d) 



IQd^l 



p ^ ad^^'^{p — a) 



2(a 



ar 



-3/2 



+ 0{i 



-5/2N 



(23) 



Hence, by using the symmetry relation (|^) and the mapping (^) onto the case of parallel 
update, after taking into account the correspondence n = i — 1 with the expansion variable 
used in [|T0|1, we obtain 



Pn = pL,r(^+ 1) 

2-p + 8d 



n 



+ 



+ 



~ 2 
(3d^/\p-(3) 



-1/2 



32v/^rfV2 4y?(/5-/3e 



n 



-3/2 + 0(^-5/2)_ 



(24) 



This result coincides with Eq. (9.22) of Ref. |jTO| except for the sign of the 0(n^^/^) term. 
The first two terms in the right-hand side of Eq. ( P^ ) generate Eq. (87) in []rT|, provided the 
slope of the density profile is considered there near the right-hand boundary of a semi-infinite 
chain, i.e., when L oo, i/L ^ 1. 



III. LOCAL DENSITY IN THE LOW-DENSITY PHASE 

The Low-density phase occupies the region a < P and a < see region A=AlUAII in 
Fig. 1. In the case of a finite chain of L sites, up to exponentially small in L corrections the 



current in the low-density (l.d.) phase is [|T2 



^^'^^-rf(a + e + r^)"p(l-a)- ^'^^ 

With respect to the asymptotic form of the density profile near the right-hand end of 
the chain, one distinguishes two subregions, AI and All. 

Subregion AI: a < and a < /3 < /3c {C, > f] > As shown in |]T2|, up to terms 



which are uniformly in z = 1, . . . , L exponentially small as L ^ oo, the local density is 
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a{l — p) 1] — T]' 



p{l 



a) 



fa + 7] + r] M 



(26) 



Hence, in the limit L ^ oo at fixed we obtain that p^\{i) is uniform, equal to the bulk 
density given by the first term in the right-hand side of Eq. (p6|). If the limit L — > cxd is 
taken at fixed j = L — i + 1, we obtain the exact expression for the non-trivial density profile 
on a semi-infinite chain with right-hand endpoint: 



„0) = ^ + 



77 — 77 



-1 



a + T] + 1] 



p{l — a) a + r] + rj ^ \a + ^ + ^ 

^/,_i(0-r//,-i(r/) 



(27) 



At large distance from the end of the chain, j ^ 1, the leading order asymptotic form 
of p^Xj) is given by the first two terms in the right-hand side of Eq. (|27|). Upon mapping 
on the case of parallel dynamics, see Eq. (|^), and taking into account the labeling corre- 
spondence j = n + 1, our result coincides with Eq. (9.25) of |jTO[. Then, the leading-order 
expression for the slope of the density profile as — > cxd reproduces exactly Eq. (82) in [11 . 

Subregion All: a < and (3 > (3c {S, > 1 and rj < 1). As shown in |T^, up to terms 
which are uniformly m i = 1, . . . , L exponentially small as L cxd, the local density in this 
subregion is 



(28) 



Hence, in the limit L — > 00 at fixed i, we obtain that p^i{i) is again uniform, equal to 
the bulk density. When the limit L — cxd is taken at fixed j = L — z + 1, we obtain the 
exact expression for the non-trivial density profile on a semi-infinite chain with right-hand 
endpoint: 



a 



p{i 



a] 



(29) 



At large distance from the right-hand endpoint, j ^ 1, the leading order asymptotic 
form of pj^r(i) subregion AH is 



pZ,Aj) 



«(i -p) 

p{l - a) 



Va + 2 
2^ 



a + 2 



3 



-3/2 



+ 0{3 



~5/2n 



(30) 



A comparison of the right-hand density profile for a finite chain, Eq. (pSj) , and its leading- 
order asymptotic form on a semi- infinite chain, Eq. (PO]), with the result of computer 
simulations for L = 300 is given in Fig. 3. 
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Upon mapping on the case of parallel dynamics, see Eq. (^), expressing the result in terms 
of the original parameters a and /?, and taking into account the labeling correspondence 
j = n + 1, we obtain 



a(l — a) a{p — a) c?^/^ 



p — a'^ 2-s/7r{p — a^) 1 — d 



ajp - a) _ Pip - 13) 



a{p — a) 



{l-a){l-dy 



n 



-3/2 



+0 



a{p — a) 



l-a)(l-d)2 



(31) 



This is the correct version of Eq. (9.28) in ||T0|; the error there occured on passing to the 
second equality in Eq. (9.27). For the slope of the density profile the first two terms in the 
right-hand side of Eq. ( |3TD yield Eq. (84) in |jTl| as a leading-order expression when n S> 1. 

The boundary between subregions AI and All: a < ac and /5 = /5c (^ > 1 and 
1] = 1). The exact expression for the density profile on a semi-infinite chain with right-hand 
endpoint is given by Eq. (p9[) at r] = 1. Its large-distance asymptotic form is redily obtained 
with the aid of expansions ([T^)- 



q(i -p) 

p{l-a) 
1 f a + 2 



-1 r'^' + o{r'^'). 



3/2 N 



(32) 



Upon mapping on the case of parallel dynamics with the aid of Eq. (^, we recover exactly 
Eq. (9.29) in dTO). Then, for the leading-order asymptotic form of the slope of the density 
profile as n ^ oo one obtains Eq. (83) in fill . 



IV. DISCUSSION 

Summarizing, we have calculated exactly the local density profiles in all the phases of the 
FASEP on semi-infinite chains with open either left-hand (particle injection) or right-hand 
(particle removal) boundary. The obtained expressions are much simpler than the finite-size 
ones, and they readily admit large-distance asymptotic analysis. 

The following general features of the local density profiles on semi-infinite chains turn 
out to be common to all the basic updates: 

(i) The profile in the maximum-current phase depends on the boundary condition at the 
finite endpoint of the chain, and does not depend on the boundary condition at infinity. The 
leading-order approach to the bulk density with the distance i from the chain end is of the 
order 0(i~^/^) and does not depend on the boundary conditions at all. 

(ii) The local density is uniform in the low-density phase of a chain with left-hand 
endpoint, and in the high-density phase of a chain with right-hand endpoint. 

(iii) The local density profiles in the low-density phase of a chain with right-hand end- 
point, and in the high-density phase of a chain with left-hand endpoint, depend on the 
boundary conditions both at the finite endpoint and at infinity. This fact is reminiscent 
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of the dependence of equilibrium states with spontaneously broken symmetry on boundary 
conditions at infinity. However, it is the continuity of the current that keeps here the in- 
formation about the boundary condition at the chain end that goes to infinity in the limit 
L — > cx). 

(iv) The analytic form of the local density profile on a semi-infinite chain with right- 
hand (left-hand) endpoint changes on passing from subregion AI (BI) to subregion All (BII) 
within the low-density (high-density) phase, compare Eqs. (^7p and (P^. The leading-order 
large-distance asymptotic approach to the bulk density is also qualitatively different: in AI 
(BI) the approach is purely exponential, with inverse correlation length = |A^^ — A~^|, 
while in All (BII) the approach is exponential, with power-law prefactor and inverse 
correlation length A^^ (A^^), where 

(33) 

On the borderline between these two subregions the leading-order asymptotic approach is 
exponential, with power-law prefactor and inverse correlation length given by A^^ 

(A~^). On crossing the mean-field line in subregion All (BII), the analytic form of the 
profile does not change, only its bending near the right-hand (left-hand) end changes from 
downward (upward) above that line to upward (downward) below it, see Eq. (^Ol). 

(v) The correlation length depends on the type of update only: for all the true discrete- 
time updates A^ is given by Eq. (|33D , and for the random sequential update, see Eq. (78) 
ini, 

We point out that from the mapping (0) of our exact results for the finite-chain current 
and local density pj^^i) 10 onto the case of parallel update, exact new representations for 



J I and pl{i) follow, which, in contrast to those found in are convenient for asymptotic 



analysis. Finally, by comparing the so obtained asymptotic results to those published in |]10 
we have corrected technical errors in two of the expressions given there. The leading-order 
large-distance expressions that follow for the slope of the density profile near the right-hand 
endpoint of a semi-infinite chain reproduce the corresponding results obtained in [ITT . 
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FIG. 1. The phase diagram in the plane of the injection and removal probabilities a and /? (see 
the text) for hopping probability p=0.75. The maximum-current phase occupies region C. Region 
A=AIUAII corresponds to the low-density phase, and region B=BIUBII to the high-density 
phase. Subrcgions Al (Bl) and All (BII) are distinguished by the different analytic form of the 
density profile. The boundary between them, P = Pc, < a < Oc {a = ac, < /3 < /3c), is shown 
by dashed segment of a straight line. The solid line a = P between subregions Al and Bl is the 
coexistence line of the low- and high-density phases. The curved dashed line is the mean-field line 
(l-a)(l -/?) = ! -p. 
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FIG. 2. The approach of the left-hand local number-density profile p versus the site nimiber 
i in the maximum-current phase {a = (3 = 0.8, p = 0.75), obtained by computer simulations for 
finite chains of length L: solid squares - L = 300, solid triangles - L = 1000 (data averaged over 
300 runs of 2^^ Monte Carlo steps per site), to the limit given by Eq. (21) - solid line. 
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FIG. 3. The behavior of the right-hand number-density profile p versus the scaled distance 
r = i/L in subregion All [a = 0.2, (3 = 0.6, p = 0.75): the solid squares present the results of 
computer simulations for a finite chain of L = 300 sites, the solid line corresponds to Eq. (29) for 
L = 300, and the dotted line - to the leading-order asymptotic form of the profile for a semi-infinite 
chain given by Eq. (30). 
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